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ABSTRACT 

Clouds, although only occupying a relatively small frac- 
tion of the tropospheric volume can have a substantial im- 
pact on the chemistry of the troposphere. In newly formed 
clouds, or in clouds with air rapidly flowing through, the 
chemistry is expected to be far more active than in aged 
clouds with stagnant air. Thus frequent cycling of air 
through shortlived clouds, i.e. cumulus clouds, is likely 
to be a much more efficient media for altering the compo- 
sition of the atmosphere than an extensive cloud cover i.e. 
frontal cloud systems. 

The impact of clouds is tested out in a 2-D channel 
model encircling the globe in a latitudinal belt form 30 
to 60° N. The model contains a detailed gas phase chem- 
istry. In addition physiochemical interactions between the 
gas and aqueous phases are included For species as H 2 0 2 , 
CH 2 0 , 0 3 and S0 2 , Henry’s law equilibria are assumed, 
whereas HN0 3 and H 2 SO 4 are regarded as completely dis- 
solved in the aqueous phase. Absorption of H0 2 and OH 
is assumed to be mass-transport limited. 

The chemistry of the aqueous phase is characterized by 
rapid cycling of odd hydrogen (H 2 0 2 , H0 2 and OH). 0 2 ~ 
(produced through dissociation of H0 2 ) reacting with dis- 
solved 0 3 is a major source of OH in the aqueous phase. 
This reaction can be a significant sink for 0 3 in the tro- 
posphere. In the interstitial cloud air, odd hydrogen is de- 
pleted, whereas NO* remains in the gas phase, thus reduc- 
ing ozone production due to the reaction between NO and 
H0 2 . Our calculations give markedly lower ozone levels 
when cloud interactions are included. This may in part ex- 
plain the overpredictions of ozone levels often experienced 
in models neglecting cloud chemical interactions. 

In the present study the existence of clouds, cloud types 
and their lifetimes are modelled as pseudo random vari- 
ables. Such pseudo random sequences are in reality dctrr- 
ministic and may, given the same starting values, be re- 
produced. The effects of cloud interactions on the overall 
chemistry of the troposphere will be discussed. In partic- 
ular, tests are performed to determine the sensitivity of 
cloud frequencies and cloud types. 


1 . INTRODUCTION. 

Measurements indicates that since the turn of the cen- 
tury there has been a marked increase in tropospheric ozone 
levels particulary over northern mid latitudes (Bojkov, 1986; 
Kley et al., 1988). The increase is a result of man made 
emissions of hydrocarbons and NO* (NO and N0 2 ) com- 
pounds (Crutzen, 1988). Tropospheric ozone is an active 
greenhouse gas, and thus increased ozone levels, particu- 
lary in the Northern Hemisphere could contribute to global 
warming. 

Several past model studies of tropospheric ozone have 
had a tendency to overpredict ozone (i.e. Isaksen et al., 
1989; Jonson and Isaksen, 1992a). Lelieveld and Crutzen 
(1990, 1991) proposed that clouds may act as an efficient 
sink for ozone, predominantly through the aqueous phase 
reaction with 0 2 . In fact, they found this reaction to 
be so fast that without a fixed ozone value, ozone levels 
would have become unrealistically low. Jonson and Isaksen 
(1992b) calculate a much smaller, but nonetheless, signif- 
icant impact of clouds on ozone levels. They found that 
about half of the ozone reductions could be attributed to 
indirect effects on the gas phase chemistry, the other half 
to the aqueous phase reaction with Oj • 

Lelieveld and Crutzen (1990, 1991) and Jonson and 
Isaksen (1992b) used crude parameterizations for clouds. 
Lelieveld and Crutzen (1990. 1991) assumed average cy- 
cling times in and out of clouds based on lagrangian statis- 
tics. Jonson and Isaksen (1992b) assumed stationary clouds 
with the air flowing through. The present study applies the 
same model, but a more flexible cloud parameterization is 
used. Based on statistics for the fractional cloud cover, 
clouds are distributed according to a pseudo random func- 
tion. 

Several species will reach concentrations in the aqueous 
phase several orders of magnitude higher than in the gas 
phase. Provided fast aqueous phase reactions occur, sev- 
eral soluble species are rapidly removed in clouds. Thus, 
important aqueous phase processes are likely to become 
far less important with increasing residence time for air in 
clouds. 
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2. DESCRIPTION. 

The model covers the latitudinal belt from 30 - 60° N. 
The horizontal resolution is 10° longitude and the vertical 
resolution 500 m from ground level to 3 km and 1 km from 3 
to 16 km, the upper boundary of the model. The timestep 
is 1 hour, and the model is run for 30 days. Mid-april 
conditions are simulated. The zonal distributions of both 
natural and anthropogenic emissions of NO* (N0-fN0 2 ), 
non-methane hydrocarbons and CO are taken from Isaksen 
et al. (1989), while emissions of S0 2 are from Jonson and 
Isaksen (1992a). For the advection, seasonally averaged 
wind fields are used (Oort, 1983). The advection equation 
is solved by an upwind scheme with corrections for numer- 
ical dispersion (Smolarkiewicz, 1983). In the vertical, the 
eddy diffusion is separated into a slow diffusive transport, 
and a rapid cloud transport. The cloud transport is a com- 
bination of direct transport from cloud base to any layer be- 
tween cloud base and cloud top, and a slower compensating 
subsidence to the layer below. This method was first de- 
scribed by Chatfield and Crutzen (1984). Cloud transport 
takes place at every longitude, regardless of whether cloudy 
or cloud free conditions, as described below, are assumed. 
However, in clouds, the cloud transport is multiplied by 
1.6, outside clouds by 0.2, empphasizing cloud advection 
where the bulk of the clouds are assumed to be concen- 
trated. The advection, eddy diffusion and cloud transport 
are described in more detail by Solberg (1989). 

Whether a grid rectangle will be in clouds (cloud mode) 
or not is determined by a pseudo random function. Each 
grid rectangle is assigned a probability for cloud based 
on statistics for the fractional cloud cover (Gordon et al., 
1984). Whenever the cloud mode is selected by the random 
function, the model has a choice of 4 cloud categories. Each 
cloud category is assigned a probability and a correspond- 
ing residence time. Probabilities for the 4 cloud categories, 
their residence times and their height intervals and optical 
thickness are listed in Table 1. If a clear sky situation is 
selected, a choice is made between 4 residence times. The 
probability for a residence time of 31 hours is 0.1, a resi- 
dence time of 6 hours, 0.2, a residence time of 2 hours, 0.4, 
and of 1 hour, 0.3. The overall probability for the indi- 
vidual cloud categories and for the length of the clear sky 
periods are determined by the product of the probabilities 
and their residence times. When selected the average resi- 
dence times for clear sky and cloud events are equeal. The 
overall probabilities for clouds are thus determined by the 
cloud fractions. This pseudo random sequence can be re- 
produced given the same starting values. Different sets of 
dissociation rates are calculated for the 4 cloud categories 
and for clear sky conditions (Jonson and Isaksen, 1992b). 

The chemistry (gas- and aqueous- phase) has been de- 
scribed by Jonson and Isaksen (1992b). In clouds, we as- 
sume that H 2 S0 4 and HN0 3 are completely absorbed by 
the droplets. For H 2 0 2) S0 2 , 0 3 , HCHO and HCOOH, 
equilibrium between the gas- and liquid- phase, determined 
by Henry’s law constants and equilibrium constants, are 


applied. The fraction of the total (gas and liquid) concen- 
tration dissolved in the aqueous phase is defined as 


' 1 + {LH e ffRT)- 1 

where L is the volume ratio between gas and liquid wa- 
ter, H e // is the efficient Henry’s constant (taking into ac- 
count the effect of dissociation), R is the universal gas 
constant and T the temperature. In clouds, the aqueous 
phase concentration is then calculated as f multiplied by 
the total (gas and liquid) concentration, and the gas phase 
concentration as (1-f) multiplied by the total concentra- 
tion. The mass-transfer of OH and H0 2 is assumed to be 
mass-transport limited, and is calculated as described by 
Schwartz (1986). In the previous study (Jonson and Isak- 





Cat. 

Prob. 

Res. 

Height (km) 

Opt. 

1 

0.1 

3 

0 - 0.5 

see caption 

2 



7 - 8(9) 

5 


0.2 

10 

5 - r 

15 




0.5(2) - 4 

20 

3 

0.2 

8 

2.5 - 5 

30 

4 

0.5 

3 

0.5(2)-6(7) 

30 


Table 1: When cloud mode is selected , Prob. denotes the 
probability for a cloud category, Res. the residence time in 
hours of the cloud event, Height , denotes the height inter - 
val(s) for the cloud layer(s) and Opt., the optical thickness 
of the cloud . For category 1 , J -values for clear sky are 
multiplied by a factor 1.5 in the lowest layer. Numbers in 
brackets denote cloud-base and cloud-top heights over land. 

sen, 1992b), we showed that aqueous phase chemistry has 
a substantial impact on the composition of the gas phase. 
The processes were studied with clouds fixed in space and 
time. This simple approach has the advantage of provid- 
ing stable conditions where the individual reactions can 
be studied in more detail. Two processes, resulting in re- 
duced ozone levels compared to gas-phase chemistry only, 
were identified: 

The first process is the aqueous phase loss through the 
reaction: 

0 3 - 0 2 + if 4 - - OH 4 20 2 (1) 

Although ozone has a low solubility in water, the reaction 
is nonetheless a significant sink for ozone. This reaction 
is also a major source of Oil in the aqueous phase. The 
sources of 0 2 ^H0 2 are absorption from the gas phase 
and recycling of hydrogen from OH through the aqueous 
phase reactions: 

OH 4 H 2 0 2 — * H0 2 - H 2 0 (2) 


1 10 


OH 4 CH 2 (OH) 2 4 0 2 — > HCOOH -f H 2 0 4 U0 2 (3) 



Long. 

fc 

short 

P 

cycling 

Qp 

Q 0 , 

fc 

long 

P 

:y cling 

Qag 

150°W 

0.63 

1.8E+6 

1.1E-6 

7.5E-7 

0.7 

2.0E+6 

1.2E-6 

6.2E-7 

80° W 

0.74 

6.8E+6 

5.0E-7 

2.3E-7 

0.69 

6.5E+6 

5.4E-7 

2.2E-7 

30°W 

0.7 

2.5E+6 

1.4E-7 

3.3E-7 

0.72 

2.4E+6 

1.5E-6 

2.4E-7 

20°E 

0.48 

6.8E+6 

5.0E-7 

2.1E-7 

0.41 

6.7E+6 

6.0E-7 

2.1E-7 

mean 

0.38 

4.3E+6 

8.7E-7 

4.4E-7 

0.39 

4.7E+6 

9.9E-7 

3.4E-7 


Table 2: Ozone production and loss terms for selected longitudes and mean values for the whole model volume, fc is 
the fractional time in the cloud mode. P is the production of ozone in molecules cm" 3 s _1 , and Q g and Q aq the gas- 
and aqueous - phase loss of ozone in s~ l averaged over a 30 day period for the lowest 6750 meters of the troposphere. 


Formic acid partially evaporates, but it also contributes to 
the cycling of hydrogen: 

OH + HCOO- + 0 2 — * C0 2 + OH- + H0 2 (4) 

OH + HCOOH + 0 2 — ♦ C0 2 + H 2 0 + H0 2 (5) 

The second process which leads to reduced ozone levels is 
due to reduced gas phase production through the reaction 
sequence: 

NO + H0 2 -+ N0 2 + OH (6) 

In sunlight N0 2 is photolyzed: 

N0 2 hu — ► NO -f~ O (7) 

O rapidly combines with 0 2 giving ozone. Whereas most of 
the odd hydrogen (OH, H0 2 and H 2 0 2 ), and thereby H0 2 , 
resides in the aqueous phase, NO* has a low solubility, and 
will remain in the gas phase. The effect is reduced ozone 
nmdiirtion. 



Figure 2: Percentage difference in ozone concentrations be- 
tween model runs without cloud chemistry and with cloud 
chemistry as a function of longitude and altitude. 


3. COMPUTATIONAL RESULTS. 



Figure 1: Diurnally averaged ozone concentrations for the 
Northern mid - latitudes as a function of latitude and alti- 
tude. Concentrations are given in parts per billion (ppb). 


In the gas phase, ozone is destroyed by a number of 
species, including H0 2 . 


( 8 ) 


In this study we are interested in quantifying the effects 
of different cloud representations rather than to study the 
individual reactions in detail. In Figure 1 ozone concen- 
trations are shown with random clouds as described above. 
Figure 2 shows the percentage difference between concen- 
trations calculated with no cloud chemistry (but with iden- 
tical dissociation rates) and calculations with cloud chem- 
istry. With cloud chemistry included, concentrations are 
lower throughout the troposphere. Ozone produced over 
polluted continents (North America, Europe and East Asia) 
is transported eastwards. During transport ozone is de- 
stroyed in clouds. In particular, large differences are found 
at the eastern rim of the Atlantic and Pacific oceans, down- 
wind from Noth America and East Asia respectively. 

Figure 3 shows the percentage difference in ozone con- 
centrations between a model run with residence times dou- 
bled, and a run where residence times have been approx- 
imately halved. Effects on ozone are generally small, and 
are slightly negative over polluted continents, most likely 
due to changes in the NO x distribution. In Table 2, ozone 
production and loss rates are shown. The tabulated values 
are averages for the whole period studied, with both cloudy 
and cloud free periods included. The effect of long and 


03 if 0 2 — ► OH 4- 20 2 






5. REFERENCES. 



LONGITUDE (DEG. W - E) 


Figure 1: Percentage difference in ozone concentrations be- 
tween model runs with long and short residence times for 
clouds. 

short cloud cycling times are compared. Production rates 
(reaction 6) and loss rates for ozone reacting with H0 2 
in both gas and aqueous phase (reaction 8 and 2 respec- 
tively) are tabulated. Over the polluted continents (80°W 
and 20° E) production rates are about a factor of 10 higher 
than over remote oceanic sites (150°W and 30° W). With 
short residence times the aqueous phase loss (reaction 2) is 
about 30% higher averaged over the whole model volume. 
This is compensated for by lower loss rates in the gas phase 
(reaction 8). However, Figure 3 show that ozone concen- 
trations generally are lower in calculations with short cloud 
cycling times. These changes are caused by indirect effects 
on the NO x chemistry and distribution rather than changes 
in aqueous phase destruction. 

4. CONCLUSIONS. 

Aqueous phase chemistry can reduce ozone concentra- 
tions in the troposphere significantly. Two mechanisms are 
identified by which ozone is reduced. Ozone absorbed in 
the droplets are consumed by a reaction with Oj . Sec- 
ondly H0 2 and NO are separated in the clouds, as H0 2 is 
absorbed by the droplets whereas NO is virtually insolu- 
ble, thus slowing ozone production. Variations in residence 
times for cloud events does not alter the ozone concentra- 
tions by more than a few percent. With such a low sensitiv- 
ity to cloud cycling parameterization, it is justified to apply 
the proposed mechanism even in models with low resolu- 
tion where individual clouds, or even entire cloud systems 
are not properly resolved. 
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